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Popular Summary: 

Dynamo models for the Sun have become increasingly more sophisticated but cannot explain the 
22-year magnetic cycle. An essential element appears to be missing. We proposed (Mayr, 

Wolff, Hartle, GRL, 28, 2001) that the periodicity of the solar magnetic cycle is determined by 
wave mean flow interactions analogous to those driving the Quasi Biennial Oscillation in the 
Earth’s atmosphere. Here we discuss some properties of the kinematic dynamo that would 
generate the magnetic fields. The toroidal magnetic field is directly driven by zonal flow and is 
relatively large in the source region Consistent with observations, this field peaks at low 
latitudes and has opposite polarities in both hemispheres. The oscillating poloidal magnetic field 
component is driven by the meridional circulation. Relative to the toroidal and the time 
independent field components, the poloidal magnetic field is small in dynamo region but is 
expected to grow as it is transported to the surface. Since the meridional and zonal flow 
oscillations are out of phase, the poloidal magnetic field peaks during times when the toroidal 
field reverses direction, which is observed. With the proposed wave driven flow oscillation, the 
magnitude of the oscillating poloidal magnetic field increases with the mean rotation rate of the 
fluid. This is consistent with the Bode-Blackett empirical scaling law, which reveals that in 
massive astrophysical bodies the magnetic moment tends to increase with the angular momentum 


of the fluid. 



Abstract: We have proposed (Mayr, Wolff, Hartle, GRL, 28, 2001) that the periodicity of the 
solar magnetic cycle is determined by wave mean flow interactions analogous to those driving 
the Quasi Biennial Oscillation in the Earth’s atmosphere. Upward propagating gravity waves 
would produce oscillating flows near the top of the radiation zone that in turn would drive a 
kinematic dynamo to generate the 22-year solar magnetic cycle. The dynamo we propose is built 
on a given time independent magnetic field B, which allows us to estimate the time dependent, 
oscillating components of the magnetic field, AB. The toroidal magnetic field, AB<p, is directly 
driven by zonal flow and is relatively large in the source region, AB</B 0 » 1 Consistent with 
observations, this field peaks at low latitudes and has opposite polarities in both hemispheres. 

The oscillating poloidal magnetic field component, ABe, is driven by the meridional circulation, 
which is difficult to assess without a numerical model that properly accounts for the solar 
atmosphere dynamics. Scale-analysis suggests that ABe is small compared to Be in the dynamo 
region. Relative to B 0 , however, the oscillating magnetic field perturbations are expected to be 
transported more rapidly upwards in the convection zone to the solar surface As a result, ABe 
(and AB«p) should grow relative to B e , so that the magnetic fields reverse at the surface as 
observed. Since the meridional and zonai flow oscillations are out of phase, the poloidal 
magnetic field peaks during times when the toroidal field reverses direction, which is observed. 
With the proposed wave driven flow oscillation, the magnitude of the oscillating poloidal 
magnetic field increases with the mean rotation rate of the fluid. This is consistent with the 
Bode-Blackett empirical scaling law, which reveals that in massive astrophysical bodies the 
magnetic moment tends to increase with the angular momentum of the fluid. 


1. Introduction 

Conventional dynamo models for the Sun have become increasingly more sophisticated. These 
models, however, suffer because they cannot explain the 22-year magnetic cycle. An essential 
element appears to be missing. 

Earlier models that placed the dynamo in the convection region failed because buoyancy 
expelled the magnetic field at a rate that was unacceptable. More recent models and helioseismic 
observations suggest that the dynamo is seated near or below the base of the convection region, 
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where the magnetic buoyancy is sufficiently small so that magnetic fields with sufficient 
magnitude can be generated (eg., Gilman, 2000). In this region of the Sun, the plasma pressure 
is so large (P » 1) that magnetic fields cannot influence the dynamics significantly. Magneto 
hydrodynamic dynamos, in which magnetic fields could counteract the flow, are not effective to 
produce the periodic 22-year oscillation. The dynamo instead must be kinematic (e.g., Chouduri, 
2000), and for that an oscillating flow is required to generate the magnetic cycle. Such an 
oscillating flow has been recently proposed by Mayr et al. (2001), and the underlying fluid 
dynamics has been elaborated by Mayr and Wolff (2001). 

The basic mechanism proposed to drive a solar flow oscillation with periodicity of 22 years 
has been studied extensively in the Earth’s atmosphere and is well understood. As was 
demonstrated in the seminal papers of Lindzen and Holton (1968) and Holton and Lindzen 
(1972), upward propagating waves, planetary waves in their case, can accelerate the zonal flow 
at low latitudes to produce a persistent periodic oscillation with a period around 24 months that 
describes the observed Quasi Biennial Oscillation (QBO). While the QBO is influenced by the 
seasonal variations of solar heating, Lindzen and Holton established that in principle such a flow 
oscillation can be generated by a steady flow of upward propagating waves without external time 
dependent forcing. 

This mechanism requires gravity waves that propagate upwards. The proposed solar dynamo 
is therefore placed in the radiative regime just below the convection region, where the dynamical 
condition of low convective stability is well suited for wave mean flow interactions. The low 
convective stability in this region helps to produce the desired oscillation period of 22 years, 
which otherwise would be much longer owing to the large length scales involved in the 
dynamics of the Sun. Employing reasonable wave amplitudes in a parameterization (Hines, 
1997a, b) that was successfully applied to describe the terrestrial QBO, a simplified analytic 
model produced oscillating zonal flow velocities of about 20 m/s at the solar equator (Mayr et 
al., 2001). With the Coriolis force vanishing at the equator, such a flow can be derived simply 
from a ID solution of the zonal momentum equation without considering the remaining 
equations for conservation of energy, mass, and momentum. In 2D modeling studies that apply a 
globally uniform wave source to the terrestrial atmosphere, the zonal flow oscillation is shown to 
be confined to low latitudes This is an important property of wave driven flow oscillations, and 
it is caused by the meridional circulation that redistributes momentum away from the equator. 
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The essential properties of this flow oscillation, such as its periodicity and the amplitudes of 
the individual flow velocities, are all determined by the wave driven fluid dynamics without 
involving the magnetic field. The flow oscillation so derived, with a 22-year cycle, would then 
drive a kinematic dynamo to generate a magnetic field oscillation with the same periodicity. In 
the following, we discuss the general properties of this dynamo. 


2. Kinematic Dynamo Generated by Flow Oscillation 

The flow oscillation to generate the dynamo would occur in the radiative interior of the Sun just 
below the convective envelop. This region is convectively stable so that gravity waves can 
propagate, and there is some observational evidence for these waves (Wolff, 1983) although they 
are difficult to detect (Palle, 1991). In Figure la we show the radial distribution of the buoyancy 
frequency, N, taken from a standard solar model. This reveals N decreasing rapidly towards the 
base of the convection region, which favors wave interactions with the background flow. Solar 
models account for radiative transfer but not dynamics, and the distribution of N is therefore 
uncertain as indicated with the dashed portion in Figure la. 

Helioseismology now provides information about the differential rotation inside the Sun. In 
Figure lb we show the best-fit solar rotation profiles of Charbonneau, et al. (1999). The 
differential rotation of the convective envelope makes the transition to the solid body rotation of 
the radiative interior through the tachocline. This reveals that the abatement of differential 
rotation extends much deeper into the interior at low latitudes than at high latitudes, although the 
rotation rate (relative to the interior) is much smaller near the equator. Fritts et al. (1998) 
proposed that waves originating in the convection region could accelerate such a flow, and it is 
understandable that it would be confined to equatorial latitudes where momentum is not 
redistributed efficiently by the meridional circulation. And in a similar vane, the wave driven 
flow oscillation we propose would also be confined to low latitudes (Mayr et al., 2001, Mayr and 
Wolff, 2001). 

In Figure lc, we illustrate what our proposed flow oscillation is like when generated by 
waves propagating out of the solar interior. Here we show a snapshot of the flow with peak 
velocity of 20 m/s and weaker reversing flow at a deeper level. The wave mechanism would 
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cause this flow pattern to migrate slowly down on a time scale of 22 years. A steady flow of 
waves would continually regenerate this reversing flow analogous to the Quasi Biennial 
Oscillation (QBO) in the terrestrial atmosphere (see Lindzen and Holton, 1968). The dashed 
curves in Figure lb illustrate how such a flow would modulate solar differential rotation during 
extremes in the solar cycle. 

In the dynamo region of the Sun, we delineate the magnetic field, B(t) = B + AB(cot), into 
time independent and time dependent components, with to = 2ji/22 years. Considering the 
electrical conductivity of the solar interior, the decay time of the magnetic field is so long that B 
may in part be the remnant of a fossil field. Whatever generates this magnetic field, we shall 
assume here that it exists and examine the time variations in AB(cot) as driven by the flow 
oscillation. 

The kinematic dynamo discussed here can be understood, to first order, in the framework of a 
2D model. We emphasize, however, that the 2D formulation for our model is a simplification 
that should be understood in the spirit of Ackhom’s razor. A 3D model would certainly be more 
realistic, but in our case, we believe, that is not essential. 

We consider the kinematic dynamo equation 

— B = V x (V x B) + q V 2 B 0) 

St 


where q is the magnetic diffusivity whose scale time is taken to be long relative to the 22-year 
cycle time in the Sun. Consequently, we neglect the corresponding term in Eq. (1). Under this 
condition, the linearized dynamo equation for the time dependent component of the magnetic 
field yields 

— AB = V x (V x AB) + V x (AV x B) ( 2 ) 

a 

with the flow, V(t) = V + AV(cot), delineated (like B) into time independent and time dependent 
components, respectively. In 2D, Eq. (2) yields the following equations for the magnetic field 
components, 
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where some of the terms, identified with brackets, will be eliminated in our simplified analysis, 
and r, 0, and cp, are respectively the radial, meridional, and zonal (azimuthal) directions. 

Given the time independent components of the magnetic field, B 9 , Be, B r , and velocities, V 9 , 
Ve, V r , the oscillating flows, AV (p (cot), AVe(cot), AV r (cot), generate in Eqs. (3) through (5) the 
oscillating magnetic field components AB^cot), ABe(cot), ABr((ot). Equation (3) describes the 
toroidal magnetic field, and Eqs. (4) and (5) describe the poloidal field, the latter being 
decoupled from the former. 

With a simplified analytic model (Mayr et al., 2001), we evaluated AV 9 (cot) at the equator 
where the Coriolis force vanishes, decoupling the zonal and meridional momentum balances. 
However, AV^oot) and AV 0 (cDt) cannot be treated without involving the full set of conservation 
equations in a 2D (or 3D) model. Notwithstanding this difficulty, we shall attempt to explore 

here some of the properties of the proposed dynamo. 

With div(V, AV) ~ 0, div(B, AB) = 0, d/dr » (l/r)5/S0, and V 9 > V e , in the oscillating flow 
region, we may assume that in general (V<p, AV,,) > (V e , AVe) > (V r> AV r ), and (B 9 , AB 9 ) > (Be, 
ABe) > (B r AB r ). Considering the known properties of the flow oscillation and resulting 
magnetic fields, we can simplify the equations by eliminating the terms in curly brackets so that 
the system becomes more tractable. It turns out that Eq (4) is then exact at the equator, where 
(B r , AB r ), (V 0 , AV 0 ) vanish. 
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3. Magnetic Fields in the Dynamo Region 


With Figure 2, we present a schematic diagram that illustrates in spherical coordinates the flow 
and magnetic field components for the proposed dynamo. 

As seen from Eq. (3), and illustrated in Figure 2, the toroidal magnetic field, AB«p, generated 
by AV,p, would peak at low latitudes where the meridional gradient of the zonal flow is largest. 
This magnetic field component would tend to peak at the height (radial distance) where the flow 
peaks. Since AB 0 would change direction at this height, we are justified in ignoring the term 
V,pABe and indicate this with a square bracket. 

In a similar vane, we can further simplify Eqs. (4) and (5). Applying Eq. (4) at the equator, 
the vertical (radial) component of the flow, AV r , interacts with the given meridional magnetic 
field, B e , to generate AB e . Since AV r tends to peak at the same height as AV„, where AB 0 
changes direction, we may ignore the term with V r AB 0 , as indicated with square bracket. 

In Eq. (5), the radial flow velocity, AV r , interacts with Be to generate AB r at mid to high latitudes. 
Since AB r and AV r are expected to peak at the same height, and AB 0 changes direction there, we 
are justified to ignore the term with V r ABe in square bracket. 

With the above simplifications, the Eqs. (3) through (5) then reduce to 

ioABm * i^AV B e 
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which describe in simplified form, and independent of each other, the toroidal magnetic field, 
AB 9 , as well as the meridional and radial components of the poloidal magnetic field, (AB 0 , AB r ). 
Equations (7) and (8) for the poloidal magnetic field satisfy div(ABe, AB r ) = 0. 

To evaluate Eqs. (6) and (7), we make the analytic ansatz that AV 9 , AV r , AB 9 vary radially 
like cos((r-r 0 ) 2nfk) in the range -X/4 < r-r 0 < X/4 as illustrated in Figure 2, where X is the vertical 
wavelength of the flow oscillation taken to be 15 Mm (Mayr et al., 2001). The corresponding 
variations for AV 0 , AB 0 , are then sin((r-r 0 ) 2ji/X), which approximately assure that there is no net 
mass transport in the meridional direction and that the magnetic flux closes (divB = 0). The 
latitude dependence of AV„ is taken to be sin((0-9O)36O/O) in the range -0/4 < 0-90 < 0/4, with 
O = 80°, so that the flow oscillation is confined within 20° around the equator. As seen from Eq. 
(6), the magnetic field, AB ?t would then have the form cos((0-9O)36O/O), centered at 20° from 
the equator, but with opposite signs in the two hemispheres. It is understood that the estimates 
we shall provide represent the amplitudes of the above analytic functions. 

With AV<p = 20 m/s, © = 2nJ22 years » 10' 8 s' 1 , r = 5xl0 8 m, and (l/r)5/50 - 4.5/r ~ 10' 8 , we 
obtain ABcp/B 0 ~ -iAV„(m/s) = -i20. The oscillating toroidal magnetic field generated by the flow 
oscillation is thus much larger than the time independent meridional component of the poloidal 
magnetic field, which assures that this field component changes direction in 1 1 year intervals. 
And the magnetic field would peak when the zonal flow changes direction (factor i). 

To estimate from Eq. (7) AB e of the poloidal magnetic field, we need to know the vertical 
(radial) flow component, AV r , which cannot be related to the wave driven zonal flow, AV<p, in a 
straight forward way. Unlike the zonal circulation at the equator, which can be derived simply 
from the zonal momentum equation, the vertical and meridional flow velocities that make up the 
meridional circulation require a 2D solution of the full set of conservation equations for energy, 
continuity and momentum. This difficulty not withstanding, the following simplified scale- 
analysis serves to provide some insight. 

Considering continuity of flow, the equations of momentum and energy conservation are 
written in the form 


iooAVtp + OAVq + — AVq, = MS(GW) (wave driven zonal momentum balance) 
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with AT the temperature perturbation around T; a the coefficient for Newtonian cooling; S the 
adiabatic temperature lapse rate or stability; Cp the specific heat at constant pressure, K the 
vertical eddy diffusivity; h = did r and L = (l/r)c?/<30 respectively the characteristic vertical and 
horizontal scales of the perturbation; Q the mean rotation rate of the solar atmosphere; and 
MS(GW) the gravity wave momentum source. For Eq. (1 1) we considered (C. L. Wolff, private 
communication) that the time constant for radiative transfer is long relative to 22 years in the 
dynamo region. 

Eliminating AT in Eqs. (10) and (1 1), we obtain the meridional and radial flow velocities 
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in terms of the azimuthal (zonal) flow velocity. 

As expected, the meridional flow, AV e , varies with the zonal flow, AV„ which is generated 
in our model by the wave interaction in Eq. (9). Substituting AV r from Eq. (12) into Eq. (7), we 
obtain ABe/B 0 = [LT/(gh 2 S)]QAV 9 . The poloidal magnetic field thus reaches peak magnitudes 
around the time when the toroidal magnetic field changes sign, which is observed at the solar 
surface. Also, the oscillating poloidal magnetic field increases with Q, consistent with the 
empirical scaling law advanced by Blackett (1947), similar to Bode s law, for massive 
astrophysical bodies that tend to reveal such a trend. 

As pointed out earlier, without a 2D numerical solution for the flow field, it is difficult to 
evaluate ABe/B e . Since AV r varies inversely with S, and the convective stability approaches zero 
at the base of the convection region as seen from Figure la, a situation develops that is 
indeterminate. Moreover, ABe/Be varies with h' 2 , and this scale parameter is not well known. In 
principle, AV r thus may take on any value, and that value could perhaps be laige enough to 
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produce in Eq. (7) a ratio ABe/B e » 1, so that the poloidal magnetic field, like the toroidal field, 
reverses even in the dynamo region. A more conservative estimate of the poloidal field would 
suggest, based on Eq. (9), that the Coriolis acceleration from the meridional flow can at most be 
comparable to the zonal flow acceleration, which means that £2AVe ~ ia>AV<p. With h = 
15x10 4 * 6 * /2ji m, L = 10 8 m, Q = 2.7X1CT 6 s' 1 , this translates for Eq. (7) into ABe/Be= -iAVe/(Lco) ~ 
A\y(LQ) = 0.1. In the dynamo region, the oscillating poloidal magnetic field component, AB e , 
is probably small compared with the time independent component. Be, and therefore much 
smaller than the toroidal magnetic field, AB r 

The above analysis, applied to the dynamo source region, may provide the boundary 
conditions for realistic 3D models that describe magnetic field transport to the solar surface. 
Upward transport of AB 9 by advection (V r term in Eq. (4) with square bracket) as well as 
convection and diffusion should be more efficient than the transport of the background magnetic 
field, B e , that presumably has a larger spatial size. In the convection region, in 3D, flux loops 
can form from the toroidal magnetic field to generate a poloidal field (a effect), as was 
demonstrated by Parker (1955, 1979) and confirmed with numerical dynamo models (see the 
reviews by Gilman, 2000; Coudhuri, 2000). One may expect then that ABe/Be, at the surface, 
increases significantly over the value we estimated for the dynamo region, so that the poloidal 
magnetic field, like the toroidal magnetic field, would reverse in 1 1-year intervals, as observed. 


4. Summary and Conclusions 

Given an oscillating flow with a period of 22 years that would develop due to wave mean flow 

interaction in the stable layer of the Sun just below the convective envelope (Mayr et al., 2001; 
Mayr and Wolff, 2001), we discussed here some properties of the kinematic dynamo that may 
generate corresponding oscillations in the magnetic field. With the kinematic dynamo being 

linear in the magnetic field, we rely on the existence of a time independent magnetic field, B, 

that would permit the flow oscillation, AV, to generate through the Lorentz force the magnetic 
field oscillation, AB. The analysis presented here was simplified in that it was restricted to the 
dynamo region. No attempt was made to explain the nature of the time independent magnetic 
field component. Also, we made no attempt to describe how the magnetic fields from the source 
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region would be transported by advection, convection and diffusion to the solar surface, which is 
of critical importance and requires a realistic 3D MHD model. 

For a given wave driven zonal flow oscillation, understood to be confined to low latitudes, 
we employed scale-analysis to estimate the associated meridional circulation and the oscillating 
toroidal and poloidal magnetic fields that are generated. Our analysis leads to the following 
conclusions: 

1) To first order, in the dynamo region, the oscillating toroidal field, AB«p, and poloidal field 
component, AB 0 , depend mainly on the time independent field component, Be. 

2) The toroidal magnetic field oscillation, AB<p, is produced primarily, and directly, by the 
oscillating zonal flow. With the wave driven zonal flow velocity, AV 9 , being largest at the 
equator, AB V peaks at low latitudes where the latitudinal gradient of the flow is largest. The 
magnetic fields in the two hemispheres have opposite polarities in agreement with observations. 
The phasing is such that AB,, peaks at the time when AV„ changes direction. Relative to B 0 , the 
oscillating component, AB,,, is large in the source region. 

3) The poloidal magnetic field oscillation is produced by the oscillating meridional 
circulation. At the equator, AB e is related to B e and to the vertical gradient of the vertical flow 
velocity. In contrast to (2), this magnetic field component is in phase with AV 9 . Since the 
meridional and zonal circulations are out of phase, the toroidal magnetic field is largest during 
times when the poloidal magnetic field changes direction, which agrees with observations. The 
ratio ABe/Be is small in the source region. 

4) Compared with the time independent component, B, the oscillating fields generated in the 
dynamo region, AB, are presumably transported more efficiently by advection, convection and 
diffusion to the solar surface. In the convection region, in 3D, flux loops can also form from the 
toroidal magnetic field to generate a poloidal field. We may expect then that the oscillations 
would dominate at the surface to produce the observed reversals in the toroidal and poloidal 
magnetic fields. 

5) Our analysis shows that wave driven flow oscillations produce a meridional circulation 
and a related poloidal magnetic field that increases with the mean rotation rate of the fluid. This 
is consistent with the empirical scaling law of Blackett (1947), also called Bode s law, which 
reveals for massive astrophysical bodies that the magnetic moment tends to increase with the 

angular momentum. 
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Figure Captions 

Figure 1 . (a) Buoyancy frequency, N, plotted versus relative radial distance to the center of the 
Sun. This is taken from a standard radiative transfer model. Since dynamics is not accounted 
for, such a model would overestimate N, and this is indicated with dashed line, (b) Azimuthal 
flow with peak velocity, U = 20 m/s, at low latitudes, which is assumed to be centered at 0.68 R 
in the tachocline and decaying flow below it. The pattern drifts slowly downward with a new 
flow developing at the top. This causes the flow at a fixed location to reverse about every 1 1 
years, (c) Over time, the solar rotation profile (solid) should shift back and forth between the 
dashed curves (assuming that the solid profile applies to the time when U = 0). This would make 
the tachocline look more, or less, extended during the 22-year cycle. (Figure taken from Mayr 
and Wolff, 2001.) 

Figure 2. Schematic that illustrates on the sphere the flow velocities, AV(o»t) = (AV 9 , AV e , AV r ) 
and the related magnetic fields, AB(cot) = (AB 9> ABe, AV r ). The toroidal field has opposite 
polarities in both hemispheres. Since the meridional and zonal (azimuthal) circulations are out 
of phase, the poloidal field peaks at the time when the toroidal field changes direction. The 
insert defines the functions that are employed in the simplified scale-analysis to describe the 
radial (height) variations of different variables. 
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